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Sondage a ponderation spatiale d'une voie de transmission 

La presente invention concerne un precede de sondage 
d'une voie de transmission formee de plusieurs canaux . II 
s'agit done d'une technique connue sous le nom de diversite 

5 de reception, technique selon laquelle un recepteur comporte 
plusieurs antennes chacune associee a ' un canal de 
transmission distinct. En d'autres termes, 1' invention 
propose une methode d' estimation de la reponse 
impulsionnelle de ces canaux de transmission. 

10 En effet, dans un systeme de transmission, notamment 

par ondes radio, le recepteur regoit pour chaque canal de 
transmission un signal emis par un emetteur. Le signal emis 
subit des fluctuations d' amplitude et de phase dans le canal 
de transmission, si bien que le signal regu par le recepteur 

15 ne lui est pas identique. Les fluctuations du signal sont 
essentiellement dues a ce que 1 ' homme de metier appelle 
1 ' interference intersymbole . Cette interference peut 
provenir de la loi de modulation employee pour la 
transmission et elle est egalement due a la propagation 

20 multi-trajets dans le canaii 

II se trouve que le signal regu est generalement issu 
d'un grand nombre de reflexions dans le canal, les 
differents trajets empruntes par le signal emis conduisant 
ainsi a des retards varies au niveau du recepteur. La 

25 reponse impulsionnelle du canal represente 1' ensemble de ces 
fluctuations, auxquelles est soumis le signal emis. II 
s'agit done la de la caracteristique fondamentale 
representant les transmissions entre 1' emetteur et le 
recepteur. 

30 La reponse impulsionnelle du canal est utilisee 

notamment par un egaliseur qui a precisement pour fonction 
de corriger 1 ' inter ference intersymbole dans le recepteur. 
Une methode classique pour realiser une estimation de cette 
reponse impulsionnelle consiste a disposer dans le signal 

35 emis une sequence d ' apprentissage formee de symboles connus . 



Cette sequence est choisie en fonction de la loi de 
modulation et de la dispersion du canal, dispersion devant 
s* entendre ici comme le retard d'un symbole emis empruntant 
le trajet le plus long du canal par rapport a ce meme 
5 symbole empruntant le trajet le plus court. La dispersion 
est couramment exprimee comme un multiple de la duree qui 
separe deux symboles emis successifs, soit un nombre de 
"duree symbole". 

A titre d'exemple, on citera deux techniques connues 
10 d' estimation de la reponse impulsionnelle d'un canal de 
transmission . 

La premiere technique fait appel a des sequences 
d' apprentissage particulieres dites sequences CAZAC, pour 
1' expression anglo-saxonne "Constant Amplitude Zero 

15 Autocorrelation". De telles sequences sont decrites dans 
1' article de A. MILEWSKI : "Periodic sequences with optimal 
properties for channel estimation and fast start-up 
equalization", IBM Journal of Research and Development, 
vol.27, N°5, Sept. 83, pages 426-431. 

20 Le systeme de radiocommunication eel lul aire numerique 

GSM fait appel a des sequences d ' apprent issage TS formees de 
2 6 symboles notes ao a a2 5 prenant la valeur +1 ou -1 . Ces 
sequences possedent les proprietes suivantes : 

20 — 
25 X = 16 

i = 5 

20 

2]a,ai^i,=0 si 0<k<5ousi 0<-k<5 

i=5 

En notant d la dispersion du canal qui vaut 4 dans le 
cas du GSM, 1' estimation de la reponse impulsionnelle prend 
30 la forme d'un vecteur X a 5 composants notes Xq a X4 . 

La sequence de symboles regrus S correspondant a la 
sequence d ' apprentissage TS est formee elle aussi de 26 
symboles notes Sq a S25. On suppose naturellement ici que 
I'emetteur et le recepteur sont parfaitement synchronises et 
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15 



dans ce cas 1 ' estimation de la reponse impulsionnelle X est 
donnee par 1' expression suivante : 



La technique CAZAC presente 1 ' avantage d'une grande 
simplicite de raise en oeuvre. Cependant, on remarque que 
chaque composante de la reponse impulsionnelle est etablie a 
partir de seulement 16 symboles regus . Etant donne que la 
sequence d ' apprentissage comprend 2 6 symboles et que la 
dispersion du canal vaut 4, il y a de 1 ' information dans le 
signal recu qui n'est pas prise en compte et cela conduit a 
une reduction des performances par rapport a 1' ideal 
theorique • 

La deuxieme technique connue fait appel au critere des 
moindres carres. Elle est decrite notamment dans les 
demandes de brevet FR 2 696 604 et EP 0 564 849. En matiere 
de rappel, cette technique fait appel a une matrice de 
mesure A construite a partir de la sequence d' apprentissage 
TS de longueur n. Cette matrice comprend (n-d) lignes et 
(d+1) colonnes, d representant toujours la dispersion du 
canal. L'element figurant a la ieme ligne et a la jieme 
colonne est le (d+i-j)ieme symbole de la sequence 
d' apprentissage : 



16 i=5 



^ 



a4 as a2~ai ao 



as a4 as azai 



a6 as a4 as az 



A = 



a 21 

J 



m — - 

4 



La sequence d ' apprentissage est choisie telle que la 
matrice A'^^A soit inversible ou 1 ' operateur . ^ represente la 
transposition. C'est naturellement le cas pour les sequences 
CAZAC mais cVest egalement le cas pour d'autres sequences. 
5 Dans la sequence de symboles regus, on ne prend pas en 

compte les quatre premiers sq a S3 car ceux-ci dependent 
egalement de symbole inconnus emis avant la sequence 
d' apprentissage, etant donne que la dispersion du canal vaut 
4. Par un abus de langage on definira done dorenavant le 

10 signal regu comme un vecteur S ayant pour composantes les 
symboles regus/ 84, S5, S5/ ... / S25* 

Des lorS/ 1' estimation de la reponse impulsionnelle 
prend la forme suivante : 

X = (At A)-l At . s 

15 Cette technique des moindres carres est un peu plus 

complexe que la precedente mais il faut noter que la matrice 
(A^ A)"^ A^ est calculee une seule fois. On remarque ici que 
chacune des composantes de 1* estimation de la reponse 
impulsionnelle X est etablie a partir de 22 symboles regus 

20 et non pas de 16 comme dans- le cas de la technique CAZAC. On 
doit done s'attendre a une amelioration des performances. 

Cependant/ quelle que soit la technique utilisee, les 
reponses impulsionnelles de chacun des canaux de la voie de 
transmission sent considerees comme independantes . 

25 La presente invention a ainsi pour premier objet un 

precede de sondage d'une voie de transmission qui prend en 
compte le fait que les differentes antennes sont 
spatialement liees . 

Selon T' invention, le procede de sondage d'une voie de 

30 transmission s' applique lorsque cette voie est formee de 
plusieurs canaux et il necessite une estimation de la 
reponse impulsionnelle C^, C^^ de ces canaux. Le 

procede comprend les etapes suivantes : 

- acquisition d'une statistique spatiale de cette voie de 
35 transmission. 



# • 
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- etablissement d'une reponse impulsionnelle corrigee au 
moins par ponderation de ces estimations de reponses 
impulsionnelles au moyen de cette statistique spatiale et 
d'une estimation du bruit additif de ces canaux. 

5 Avantageusement, cette statistique spatiale correspond 

a une estimation de la correlation des canaux de 
transmission pris deux a deux. 

Suivant un mode de realisation pref erentiel , cette 
estimation de la correlation des canaux de transmission 

10 prend la forme d'une matrice de correlation spatiale ainsi 
d^finie que 1' element figurant a la ieme ligne et a la jieme 
colonne soit obtenu par lissage du produit de la transposee 
hermitienne de la reponse impulsionnelle estimee du ieme 
canal et de la reponse impulsionnelle estimee du jieme 

15 canal. 

Selon une caracteristique additionnelle, un signal 
recu S par un premier canal correspondant a une sequence 
d'apprentissage emise, en notant A la matrice de mesure 
associee a cette sequence d'apprentissage, 1' estimation du 
20 bruit additif Nqi de ce canal est obtenue par normalisation 
de I'energie du vecteur (S - ACi). 

On peut prevoir que cette normalisation soit suivie 

d'une.etape de moyennage. 

De plus, une matrice de bruit N etant formee k partir 
25 de 1' estimation du bruitT additif Nqi, No2/ ^On 

canaux, une matrice de ponderation spatiale G' etant ainsi 
definie a partir de la matrice de correlation spatiale G et 
de la matrice de bruit : 

S- = G (G + N)-l, 
30 les reponses impulsionnelles corrigees C'l, C'2f / C'n 

sont obtenues par 1' expression suivante : 
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Le precede de sondage d'une voie de transmission est 
done mis en oeuvre a partir de 1 ' estimation de la reponse 
impulsionnelle des differents canaux consideres comme autant 
5 de canaux independants . Or les erreurs d' estimation sont 
inevitables* De maniere generale, la determination de la 
reponse impulsionnelle d'un canal de transmission unique est 
un probleme qui ne peut etre resolu de fagon exacte en 
presence de bruit additif. De plus les techniques 

10 anterieures font 1 ' hypothese implicite que cette reponse 
impulsionnelle peut prendre une forme queleonque. 

La pre sent e invention a ainsi pour deuxieme objet une 
methode d' estimation de la reponse impulsionnelle d'un canal 
de transmission qui presente une meilleure resistance au 

15 bruit additif ou, autrement dit, qui mene a une erreur 
reduite comparee a 1 ' erreur d' estimation des techniques 
connues. Cette methode sera avantageusement appliquee pour 
I'un au moins des canaux formant la voie de transmission 
avant 1 ' etablissement des reponses impulsionnelles corrigees 

20 par ponderation au moyen de la statistique spatiale et d'une 
estimation du bruit additif de ces canaux. 

Selon 1' invention, cette methode necessite un signal 
recu par un canal, ce signal regu correspondant a une 
sequence d' apprentissage emise. La methode comprend les 

25 e tapes suivantes : 

- acquisition d' une statistique temporelle de ce canal de 
transmission, 

- etablissement de 1' estimation de la reponse impulsionnelle 
de ce canal, cette estimation etant ponderee par la 

30 statistique temporelle du canal au moyen du signal regu. 
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La statistique teinporelle du canal represents une 
valeur de la reponse impulsionnelle anterieure a 
1* acquisition du signal regu. La susdite ponderation 
introduit le fait que la reponse impulsionnelle afferente au 
5 signal regu a une valeur probablement plus proche de cette 
valeur anterieure qu'une valeur qui en serait tres eloignee. 
Ainsi, statistiquement, I'erreur d' estimation est diminuee. 

Avantageusement, cette statistique correspond a une 
estimation de la covariance de cette reponse impulsionnelle • 
10 Selon une premiere variante du procede-, celui-ci 

comprend les etapes suivantes : 

- lissage de la reponse impulsionnelle et orthonormalisation 
au moyen d'une matrice de transformation W pour obtenir 
1' estimation de la covariance qui prend alors la forme d'une 

15 matrice L' , 

- recherche des vecteurs propres Vi ' et valeurs propres A-i' 
associees de cette matrice L', 

- estimation de la reponse impulsionnelle instantanee du 
canal a partir du signal regu et application de cette 

20 matrice de transformation W.. pour former un vecteur X', 

1 ' etablissement de 1' estimation ponderee Xp etant ainsi 
realisee : 



ou Nq est un nombre reel strictement positif representant le 
25 bruit additif du canal. 

On peut ici prevoir que le bruit additif soit choisi 
egal a la plus petite des valeurs propres Xj_ ' . 

Lorsqu^'un sous-ensemble de ces valeurs propres Xi ' 
presente une contribution inferieure a un seuil 
30 predetermine, on peut egalement prevoir que chacune de ces 
valeurs propres soit forcee a la valeur du bruit additif. 

La complexite s ' en trouve reduite d'autant. 

Selon une seconde variante du procede, 1' estimation de 
la covariance prenant la forme d'une matrice R, en notant A 
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la matrice de mesure associee- a la sequence d • apprentissage, 
I'etablissement de 1' estimation ponderee est ainsi realise : 

Xp = (A^A + NqR"!)-! At.s 
ou No est un reel strictement positif representant le bruit 
5 additif du canal. 

Par ailleurs, il est possible d'effectuer une 4tape 
d'orthonormalisation de la matrice R au moyen d'une matrice 
de transformation W pour obtenir une nouvelle matrice R', 
1' estimation ponderee prenant alors la nouvelle forml 
0 suivante : 

Xp = W(I + NqR'"!)-! A't.s 
ou la matrice A' est egale au produit de la matrice de 
transformation W et de la matrice de mesure A. 

Avantageusement, 1' expression (i + NqR'-I)"! est 
5 calculee au moyen du lemme d'inversion matricielle. 

La presente invention apparaitra maintenant de maniere 
plus detainee dans le cadre de la description qui suit ou 
sont proposes des exemples de mise en oeuvre a titre 
illustratif, ceci en reference aux figures annexees qui 
0 representent : 

la figure 1, un diagramme identifiant les 
principales etapes d'un mode de mise en oeuvre du precede 
de sondage d'une voie de transmission selon 1 ' invention, 

- la figure 2, une jsremiere variante d' estimation de 
la reponse impulsionnelle d'un canal selon 1' invention, 

la figure 3, une deuxieme variante d' estimation de 
la reponse impulsionnelle d'un canal selon 1' invention. 

L' invention sera presentee dans son application au GSM 
car ce systeme a le merite d'etre bien connu de I'homme du 
metier. II s'agit done la d'une presentation adoptee dans un 
souci de clarte mais il ne faut y voir en aucun cas une 
limitation de 1' invention a ce seul systeme. 

En reference a la figure 1, le procede de sondage 
d'une voie de transmission s' applique lorsque cette voie 
comporte au moins deux canaux de transmission, n canaux dans 
le cas general. Les canaux correspondent chacun a une 
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antenne distincte. On considere done que chaque canal a fait 
I'objet d'une estimation de sa reponse impulsionnelle 
respective Ci, Cn selon 1 ' une quelconque des 

techniques disponibles ♦ 

En premier lieu le procede prevoit 1 ' acquisition d'une 
statistique spatiale de cette voie de transmission. Par 
statistique spatiale on entend un ensemble de donnees 
refletant le comportement de cette voie sur une periode 
predeterminee que I'on conviendra d'appeler periode de 
correlation. En effet, les differentes antennes etant fixes, 
les signaux regus sur ces antennes presentent un certain 
degre de correlation. L' invention entend precisement tirer 
parti de ce fait pour ameliorer la qualite d' estimation de 
la reponse impulsionnelle de 1 ' un des canaux au moins. A 
titre d'exemple, une telle statistique peut etre obtenue au 
moyen d'une matrice de correlation spatiale G : 

v(c!;ci) v(cjc2) v(c!;c„) ^ 



G = 



V(c5Ci) V(C5C2) V(C5C„) 



U(c:;co v(c|;c2) v(c:;c,) ; 

ou I'operateur represente la transposition hermitienne. 

Cette matrice G carree de dimension (n,n) peut ainsi 
etre representee de maniere generique par 1' element g^j 
figurant a la ieme ligne et a la jieme colonne : 
g-. = V(C^Cj) 

Get element gij est obtenu par lissage du produit C^Cj 
au moyen des reponses impulsionnelles estimees Ci, Cj des 
ieme et jieme canaux obtenues durant la periode de 
correlation. Ce lissage est une estimation de la correlation 
de ces deux canaux. 

On entend ici lissage dans un sens tres general, 
c'est-a-dire toute operation permettant de lisser ou de 
moyenner le produit C^Cj sur la periode de correlation. 
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Un premier exemple de lissage consiste a faire la 
moyenne de ce produit sur cette periode supposee comprendre 
e estimations successives de chacune des reponses 
impulsionnelles Ci, C2, Cn : 



Un deuxieme exemple de lissage consiste ^ actualiser, 
a la pieme estimation recue pour chacun des ieme et jieme 
canaux, I'expression de lissage Vp_i(C^Cj) obtenue a la 
10 (p - l)ieme estimation au moyen d'un coefficient 
multiplicatif a, ce coefficient compris entre 0 et 1 etant 
souvent connu sous le nom de facteur d'oubli de lissage : 
Vp(c^Cj) = aC?Cj. + {l-a)Vp.,(C?Cj) 

L' initialisation peut se faire par tous moyens, 
15 notamment au moyen de la premiere estimation obtenue ou bien 
par une moyenne obtenue comma dans le premier exemple pour 
les premieres estimations reques . 

Le precede de sondage propose ensuite de realiser une 
estimation du bruit additif Nqi, N02, Non present dans 

!0 chacun des canaux au moyen des estimations des reponses 
impulsionnelles respectives Ci, C2, de ces canaux. 

II sera propose diff§rentes solutions pour mener A 
bien cette estimation du bruit et on traitera le cas d'un 
seul canal, le premier phr exemple, etant donne que les 
:5 memes solutions s'appliquent a chacun des canaux. 

Une solution simple consiste a affecter Nqi d'une 
valeur predeterminee qui reflete un seuil en dessous duquel 
il est peu probable que le bruit additif puisse descendre. 
Cette valeur pourrait etre determinee par une mesure de 
0 rapport signal a bruit, ou par les performances du 
recepteur, ceci a titre d' exemple. 

Par ailleurs, 1' estimation du bruit Nqi peut etre 
obtenue au moyen de 1' estimation de la reponse 
impulsionnelle Ci du premier canal et du signal regu 
5 correspondant S. Cette estimation du bruit a done pu etre 
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realisee avant la mise en oeuvre du precede de sondage a 
ponderation spatiale selon quelque technique que ce soit. 
Toutefois, si ce n'etait pas le cas, on propose ici une 
methode appropriee lorsque 1' estimation de la reponse 
5 iitipulsionnelle Ci est acquise selon la technique des 
moindres carres. 

II s'agit done d'estimer directement le bruit additif 
a partir du signal regu S et de la inatrice de mesure A, En 
effet/ en notant le vecteur bruit affectant le signal 

10 regu, il vient que : 

S = AC]_ + Ni 

Compte tenu du fait que les vecteurs S et Ni ont 22 
composantes, 1 ' estimation du bruit Nqi prend la forme 
suivante .: 

15 Noi - (£^) (S - ACi)h (S - ACi) 

Naturellement cette estimation du bruit additif Nqi 
peut etre moyennee ou lissee. 

A partir des estimations du bruit additif des 
differents canaux Nqi. No2/ • — / ^On^ ^n construit une 
20 matrice de bruit N dont 1' element figurant a la ieme ligne 
et a la jieme colonne est note n^ j . 

Lorsque i est egal a j, cet element nij prend la 
valeur du bruit estime Nqi. 

Lorsque i et j sont -differents, 1' element n^j est nul . 
25 Par ailleurs, toujours dans un souci de 

simplification, il est possible de calculer la valeur 
moyenne Nq des bruits estimes Nqi/ ^02/ ^On de 

forcer chacun des elements n^j de la diagonale de cette 
matrice N a cette valeur moyenne Nq . En notant I la matrice 
30 identite, il s'ensuit que la matrice de bruit N prend la 
forme suivante : 

N = Nq I 

La matrice de correlation G et la matrice de bruit N 
permettent de definir une nouvelle matrice, la matrice de 
35 ponderation spatiale G' : 
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G' = G (G + N) -1- 
Les reponses impulsionnelles corrigees C i, C'2/ 
C'n moyen de cette matrice de ponderation spatiale G' 

sont ainsi definies : 



^ 2 



C • * / 



^2 



= G 



ou 1' operateur .t represente toujours la transposition. 

L'une au moins des reponses impulsionnelles corrigees 
C'l est utilisee a la place de 1' estimation de la reponse 
impulsionnelle dans le recepteur. 

A titre optionnel, 1' invention permet egalement 
d'ameliorer 1' estimation de la reponse impulsionnelle C^, 
C2/ ..•/ Cn de chacun des canaux, ceci avant la mise en 
oeuvre du precede de sondage a ponderation spatiale d'une 
voie de transmission, c*est-a-dire avant 1 ' etablissement de 
la ou des reponses impulsionnelles corrigees C'l, 



n ' 



On prevoit ainsi 1 ' acquisition d'une statistique 
temporelle pour 1 ' un au moins de ces canaux, le premier par 
exemple. Par statistique temporelle, on entend un ensemble 
de donnees refletant le comportement de ce canal considere 
independamment des autres canaux sur une periode 
predeterminee que 1 ' on conviendra d'appeler periode 
d' analyse. II s'agit done d'une representation du 
comportement moyen du canal pendant la periode d' analyse. 
Cette statistique peut etre etablie par quelque moyen que ce 
soit et en quelque lieu que ce soit. En effet, 
I'etablissement de la statistique peut prendre place dans un 
equipement quelconque du reseau de radiocommunication, Ce 
qui importe c'est que le recepteur puisse acquerir cette 
statistique. 
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A titre d' example, une telle statistique temporelle 
peut etre obtenue de la maniere suivante . 

Au moyen du signal regu S pendant la periode 
d' analyse, on calcule une estimation X de la reponse 

5 impulsionnelle selon une methode connue. 

Si I'on retient la technique des moindres carres cette 
estimation X vaut : 

X = (At A)-l At.s (1) 
On rappelle ici que I'emetteur et le recepteur sont 

10 supposes synchronises a mieux que un demi-symbole pres, 
auquel cas le signal de reception est forme par le vecteur S 
dont les composantes sont les symboles regus S4 a S25 
synchrones des symboles a4 a a25 1^ sequence 

d' apprentissage TS . Si une telle synchronisation n'etait pas 

15 acquise, plusieurs solutions sont disponibles pour 
I'acquerir et on en citera deux a titre d'exemple. 

La premiere solution consiste a decaler le signal regu 
en avance ou en retard de j periodes symboles de sorte que 
Sjt =: (s4-j, S5-j, sg-j, S25-j)- 

20 On calcule alors 1' estimation Xj pour chacun des 

vecteurs Sj et on retient la valeur pour laquelle Xj^-Xj 
est maximal. Cette valeur donne la synchronisation 

escomptee et il suffit de remplacer le vecteur S dans 
1' equation (1) par le vecteur SjM- 

25 La deuxieme solution consiste a augmenter 

artif iciellement la dispersion d du canal d'une quantite 2q 
predeterminee . On peut alors definir une matrice de mesure 
modifiee Aj^ comprenant (n-d-2q) lignes et (d+2q+l) colonnes. 
En reprenant' les valeurs 26 et 4 respectivement pour n et 

30 d : 
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A m ~ 



V 



a 4+2q ••••a4 a 33231 ao 
a 5+2q '•••a5a4a3a2ai 
a 6+2q ••'•a6a5a4a3a2 

a 7^2q 

a 25 a 21-2q 
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15 



20 



25 



II faut alors 

reduire le nombre de composantes du signal regu S de cette 
meme quantite 2q et on peut convenir de retenir le vecteur 
Sj^ modi fie : 

(s'4+q, s'5+q, . . . s'25-q) 



Sm^ = 



On obtient ainsi une estimation modifiee X 



>^m = (Ain^ Ain) K^'S^ 
Cette estimation modifiee 
composantes : 



m 



^m^ ~ (^-a/ 



^O' ^1^ 



(jn comporte d+2q+l 



X4 , ... X4^q) 



On recherche alors la valeur jj^ de j comprise entre -q 
et +q qui maximise 1' expression suivante : 

2 X J + k X j+k 

k = 0 — 

ou I'operateur .* represente la conjugaison complexe, 

Cette valeur jj^ determine 1' estimation X de la reponse 
impulsionnelle pour une dispersion d=4 : 

en 



La synchronisation s ' en deduit immediatement 
appliquant le decalage jj^ au signal regu S. 

On peut des lors appliquer a nouveau I'equation (1). 

On construit maintenant une matrice de lissage L par 
lissage des differentes estimations X obtenues pendant la 
periode d* analyse pour obtenir une estimation de la 
covariance associee a cette reponse impulsionnelle. On 
entend ici lissage dans un sens tres general, c'est-a-dire 
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toute operation permettant de lisser ou de moyenner la 
reponse impulsionnelle sur la periode d' analyse. On obtient 
ainsi une representation statistique du comportement du 
canal de transmission. Ce lissage peut etre obtenu notamment 
5 par I'une des deux methodes proposees ci-dessus, 
1' expression la plus simple de cette matrice de lissage L 
etant la suivante : 

ou m correspond au nombre de sequences d' apprentissage sur 
10 lequel est calcule le lissage. 

On admet ici que cette matrice de lissage peut etre 
approchee par 1' equation suivante : 

L r (A^A) "1 No + R (2) 
ou Nq represente toujours le bruit present dans le 
15 canal de transmission ou bruit additif et ou R est une 
matrice que 1 ' on a coutume d'appeler statistique a priori du 
canal car elle represente le comportement du canal 
abstraction faite du bruit. 

On admet egalement que la matrice de mesure A est bien 
20 conditionnee, c'est-a-dire que les valeurs propres de la 
matrice A^a sont tres proches les unes des autres. Dans ce 
cas, il est interessant de proceder a 1 ' orthonormalisation 
des vecteurs constitues par les lignes de la matrice de 
mesure A, mais il ne faut pas voir la une limitation de 
25 1' invention. 

Pour ce faire, on emploie une matrice de 
transformation W telle que : 

A'- = AW et A'^A' = I 
oil I represente la matrice identite. 
30 En notant L' la matrice ainsi definie : 

L = WL'Wt, 

il vient que Inequation (2) peut maintenant s'ecrire : 
L' - NqI + R' (3) 
Selon une premiere variante representee a la figure 2, 
35 on remarque que les vecteurs propres vi ' de L* et Vi de R' 
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z identiques tandis que les valeurs propres ' de L' et 
de R' sont decalees de Nq . Soit en prenant toujours la 
e valeur de 4 pour la dispersion du canal, pour tout i 
cris entre 0 et 4 : 

Vi' = Vi 

Xi' = Xi + No 

II apparait ainsi que la determination des vecteurs 
:res et valeurs propres de R' et celle de L' sont 
""iques sous reserve que Nq soit ccnnu. 

L'etape d'estimation du bruit sera decrite plus loin 
r rendre 1' expose plus clair, bien que cette etape 
=cede celle qui va maintenant etre explicitee. 

II convient done maintenant de rechercher les couples 
:ur-propre/vecteur-propre pour 1 ' une ou 1 ' autre des 
rices L' ou R' . Cette etape ne sera pas plus detaillee 
bien connue de I'homme du metier. Par ailleurs, il va 
: dire que 1 ' on peut annuler les valeurs propres dont la 
:ribution est jugee non significative. Par exemple, si 
valeurs propres sont classees par ordre ' decroissant , on 
:rime les dernieres qui- sont telles que leur somme soit 
rrieure a un seuil predetermine. 

L' etape suivante consiste a estimer la reponse 
,^lsionnelle instantanee X etablie selon 1 ' une quelconque 
techniques connues a partir du signal regu correspondant 
a derniere sequence "ci' apprentissage regue. En notant 
WX \ cette derniere estimation est ponderee par la 
.ode suivante pour obtenir une ponderation temporelle Xp 
a reponse impulsionnelle instantanee : 



.=0 ^ ^0 J 



4 r 



Pour obtenir la ponderation Xp il faut done estimer le 
t addi tif Nq . 

Cette estimation du bruit peut etre effectuee 
minent selon 1 ' une des methodes exposees plus haut. 



Une autre possibilite consiste a considerer que la 
derniere valeur propre, (la plus faible) de la matrice de 
lissage L est egale a Nq ; 

X^' = Nq ou = 0. 
5 Quelle que soit la methode retenue/ la ponderation 

temporelle Xp de 1' estimation de la reponse impulsionnelle 
instantanee peut alors etre realisee comme mentionne ci- 
dessus . 

Selon une seconde variante representee a la figure 3, 

10 on etablit 1 ' estimation ponderee Xp directement comme suit : 

Xp = (A^A + NoR"^)"^At.S 

ou bien en reprenant la matrice de transformation W 

def inie ci-dessus : 

Xp = W (I + NqR'"^)"^ A't.s (4) 

15 Conformement a 1' equation (3) : 

R' = L' - NqI 

II convient la aussi d'estimer le bruit additif Nq. 

Une solution avantageuse pour obtenir la ponderation 

temporelle Xp consiste a adopter la methode suivante . 

20 On divise la matrice R' par Nq : 

R ' 

B = — 

N 0 

II s'ensuit que : 

I -I- NqR'"^ = I + B~l 
On utilise le lemme d' inversion matricielle pour 
25 calculer la matrice de ponderation P = (I + B"*l)~^. 

Ainsi, en notant e^ les vecteurs canoniques, on 
procede a 1' iteration suivante : 
- initialisation : 
P = B 

30 - pour i variant de 0 a d (4 dans le cas present) : 

Pe, (Pe^ )^ 



P = P - 



1 + e^Pe^ 



P etant connu, il reste a etablir la ponderation Xp 
selon Inequation (4). 
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On remarquera que la matrice de ponderation P n'est 
pas necessairement calculee a 1' occasion de chaque nouvelle 
sequence d' apprentissage emise . Elle peut etre calculee a un 
rythme plus lent car elle varie sensiblement au meme rythme 
5 que R' et done plus lenteinent que le signal reQU S, 

On remarquera egalement que 1 ' estimation ponderee est 
realisee sans avoir recours a la reponse impulsionnelle 
instantanee. Elle est produite directement a partir du 
signal reQU S. 

10 Quelle que soit la variante retenue, c'est done 

1' estimation ponderee Xp qui sera avantageusement utilisee 
comitie estimation de la reponse impulsionnelle pour la 

mise en oeuvre du precede de sondage a ponderation spatiale 
d'une voie de transmission, c'est-a-dire pour produire une 

15 ou plusieurs reponses impulsionnelles corrigees. 
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REVEND I CATIONS 

1) Procede de sondage d'une voie de transmission 
formee de plusieurs canaux, ce procede necessitant une 
estimation de la reponse impulsionnelle (Ci, C2/ /^n) de 

5 ces canaux, caracterise en ce qu'il comprend les etapes 
suivantes : 

- acquisition d'une statistique spatiale de cette voie de 
transmission, 

- etablissement d'une reponse impulsionnelle corrigee (C'l, 
10 C'2' '^'n) moins par ponderation desdites estimations 

de reponses impulsionnelles (Ci, C2, /Cn) au moyen de 

ladite statistique spatiale et d'une estimation du bruit 
additif (Nqi/ No2/ ^On) desdits canaux. 

2) Procede selon la revendication 1, caracterise en ce 
15 que ladite statistique spatiale correspond a une estimation 

de la correlation desdits canaux de transmission pris deux a 
deux. 

3) Procede selon la revendication 2, caracterise en 
ce que cette estimation de la correlation des canaux de 

20 transmission prend la forme d'une mat rice de correlation 
spatiale (G) ainsi definie que 1' element (gij) figurant a la 
ieme ligne et a la jieme colonne soit obtenu par lissage du 
produit (C-Cj) de la transposee hermit ienne de la reponse 

impulsionnelle estimee (Ci) du ieme canal et de la reponse 
25 impulsionnelle estimee (Cj~) du jieme canal. 

4) Procede selon la revendication 3 caracterise en 
ce que, un signal regu S par un canal correspondant a une 
sequence d' apprent issage emise, en notant A la matrice de 
mesure assdciee a ladite sequence d' apprentissage/ 

30 1' estimation du bruit additif (Nqi) de ce canal est obtenue 
par normalisation de I'energie du vecteur (S - ACi) . 

5) Procede selon la revendication 4, caracterise en ce 
que ladite normalisation est suivie d'une etape de 
moyennage . 

35 6) Procede selon I'une quelconque des 

revendications 3 a 5 caracterise en ce que, une matrice de 
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bruit (N) etant formee a partir de 1 ' estimation du bruit 
additif (Nqi/ No2' ^on) canaux, une matrice de 

ponderation spatiale (G') etant ainsi definie a partir de 
ladite matrice de correlation spatiale (G) et de ladite 
5 matrice de bruit : 

G' = G (G + N)"l, 
lesdites reponses impulsionnelles corrigees (C'x' C*2' ' 
C'n) sent obtenues par 1' expression suivante : 





1 t 

1 


\ 








c 


1 t 

2 
















= G' 
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10 7) Precede selon I'une quelconque des 

revendications 1 a 6 caracterise en ce que, avant 
1 ' etablissement de ladite reponse impulsionnelle corrigee 
(C'l) d'un canal, le signal regu (S) par ce canal 
correspondant a une sequence d ' apprentissage emise, il 

15 comprend les etapes suivantes : 

- acquisition d'une statistique temporelle de ce canal de 
transmission, 

- etablissement de 1' estimation de la reponse impulsionnelle 
(Ci) dudit canal, cette estimation (Xp) etant ponderee par 

20 la dite statistique temporelle du canal au moyen dudit 
signal regu (S) . 

8) Precede selon la revendication 7, caracterise en 
ce que ladite statistique temporelle correspond a une 
estimation de la covariance de ladite reponse 

25 impulsionnelle . 

9) Procede selon la revendication 8 caracterise en 
ce qu*il comprend les etapes suivantes : 

lissage de ladite reponse impulsionnelle et 
orthonormalisation au moyen d'une matrice de transformation 
30 W pour obtenir ladite estimation de la covariance qui prend 
alors la forme d'une matrice L', 
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- recherche des vecteurs propres (vi') et valeurs propres 
(Xi') associees de cette matrice L' , 

- estimation de la reponse impulsionnelle instantanee du 
canal a partir dudit signal recu (S) et application de cette 

5 matrice de transformation W pour former un vecteur X', 

1 'etablissement de ladite estimation ponderee (Xp) etant 
ainsi realisee : 

Xp = Z[^^nvtX')]wv' 

oti Nq est un nombre reel strictement positif representant le 
10 bruit additif dudit canal. 

10) Precede selon la revendication 9, caracterise en 
ce que ledit bruit additif (Nq) est choisi egal a la plus 
petite desdites valeurs propres {^i')- 

11) Precede selon 1 ' une quelconque des revendications 9 
15 ou 10, caracterise en ce qu'un sous-ensemble desdites 

valeurs propres (^i') presentant une contribution inferieure 
a un seuil predetermine, chacune de ces valeurs propres est 
forcee a la valeur dudit bruit additif (Nq) . 

12) Precede selon la revendication 8 caracterise en 
20 ce que ladite estimation de la covariance prenant la forme 

d'une matrice R, en notant A la matrice de mesure associee a 
ladite sequence d ' apprentissage, 1 ' etablissement de ladite 
estimation ponderee (Xp) est ainsi realise : 
Xp = (A^A + NqR"^)"^ A^.S 
25 ou Nq est un reel strictement positif representant le bruit 
additif dudit canal. 

13) Precede selon la revendication 12, caracterise 
en ce qu'il" comprend une etape d' orthenormalisation de 
ladite matrice R au meyen d'une matrice de transformation W 

30 pour obtenir une nouvelle matrice R', cette estimation 
ponderee prenant alors la nouvelle forme suivante : 

Xp = W(I + NqR'"^)"^ A*".S 
ou la matrice A' est egale au produit de cette matrice de 
transformation W et de ladite matrice de mesure A. 
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14) Precede selon la revendication 13, caracterise 
en ce que 1' expression (I+NqR'"^)"^ est calculee au moyen du 
lemicie d' inversion matricielle. 
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